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Incremental  Length  Diffraction 
Coefficients  for  the  Shadow  Boundary 

of  a  General  Cylinder 


1  Introduction 

In  the  physical  theory  of  diffraction  (PTD)  [1],  the  scattered  field  is  divided  into  a  physical 
optics  (PO)  field  and  a  nonuniform  (NU)  field.  The  PO  field  for  perfectly  electrically  conducting 
(PEC)  objects  is  obtained  from  an  integration  of  the  PO  current  over  the  object.  One  reason  for 
the  limited  accuracy  of  the  PO  field  is  that  the  PO  current  fails  to  closely  approximate  the  exact 
current  in  the  vicinity  of  shadow  boundaries  and  on  the  shadow  side.  A  significant  improvement 
in  the  accuracy  of  the  computed  fields  can  therefore  be  obtained  by  finding  good  approximations 
for  the  fields  radiated  by  the  NU  currents  near  the  shadow  boundary.  Approximations  for  these 
NU  fields  can  be  obtained  by  first  integrating  the  product  of  the  free-space  Green’s  function  and 
the  approximated  NU  currents  excited  on  a  strip  of  differential  width  transverse  to  the  shadow 
boundary  of  a  canonical  two-dimensional  (2D)  scatterer  that  closely  conforms  locally  to  the  shape 
of  the  actual  scatterer  in  the  vicinity  of  the  shadow  boundary.  The  differential  fields  obtained 
by  this  integration  are  known  as  the  incremental  length  diffraction  coefficients  (ILDC’s).  Once 
the  ILDC’s  are  determined,  the  NU  fields  are  obtained  by  integrating  them  along  the  shadow 
boundary  of  the  actual  scatterer. 

Although  for  some  time  ILDC’s  have  been  available  for  the  shadow  boundaries  of  two- 
dimensional  (2D)  PEC  scatterers  with  sharp  edges  such  as  the  wedge  [2],  [3],  [4],  [5]  it  is  only 
recently  that  they  have  been  obtained  for  convex,  smoothly-curved  cylinders.  Hansen  and  Shore 
obtained  ILDC’s  for  the  shadow  boundary  of  a  circular  cylinder  [6]  and  a  parabolic  cylinder  [7] 
using  an  earlier  form  of  the  method  of  this  report,  and  presented  the  basic  ideas  of  this  report  in 
[8].  Yaghjian  et  al.  [9]  obtained  shadow-boundary  ILDC’s  for  a  convex  cylinder  by  approximat¬ 
ing  the  far  fields  radiated  by  the  NU  shadow-boundary  currents  of  a  convex  cylinder  by  the  far 
fields  radiated  by  the  NU  shadow-boundary  currents  of  circular  cylinders,  and  then  substituting 
the  approximated  fields  into  general  expressions  [5],  [9],  [10].  These  approximate  ILDC’s  were 

received  for  publication  1  Aug  1997 


1 


then  modified  to  account  for  a  varying  radius  of  curvature  in  the  shadow  region. 

The  purpose  of  this  report  is  to  show  how  ILDC’s  can  be  obtained  for  the  NU  currents  near 
the  shadow  boundaries  of  PEC  2D  convex  cylinders  of  general  cross-section  and  smoothly  varying 
radius  of  curvature,  by  approximating  and  integrating  the  NU  currents.  We  derive  asymptotic 
expressions  for  the  ILDC’s  that  are  rapidly  calculable  and  usable  in  general-purpose  computer 
codes.  This  method  of  obtaining  convex-cylinder  shadow-boundary  ILDC’s  is  an  alternative  to 
the  field  substitution  of  Yaghjian  et  al.  [9]  described  above.  The  field  substitution  method, 
because  it  works  directly  with  fields,  is  perhaps  simpler  than  the  method  of  this  report.  For 
applications  in  which  heightened  accuracy  is  important,  however,  the  NU  current  approximation 
and  integration  method  described  here  may  be  preferable  because  it  yields  more  accurate  ILDC’s 
than  does  the  field  substitution  method  [11  j. 

The  report  is  organized  as  follows.  In  Section  2  we  obtain  the  Fock  approximations  of  the 
NU  currents  on  a  convex  cylinder.  These  approximations  are  validated  by  specializing  them 
to  the  cases  of  a  circular  and  parabolic  cylinder  and  showing  that  the  Fock  approximations  to 
the  NU  currents  closely  match  the  exact  NU  currents.  In  Section  3  we  show  how  the  ILDC’s 
corresponding  to  the  Fock  current  approximations  of  Section  2  can  be  obtained  efficiently  by 
using  an  accurate  and  rapidly  calculable  approximation  to  the  integral  defining  the  ILDC.  In 
Section  4  we  show  how  to  determine  the  canonical  2D  cylinders  from  the  actual  3D  scatterers. 
As  an  example  of  obtaining  and  using  shadow  boundary  ILDC’s  we  show  results  obtained  for 
the  field  scattered  by  a  PEC  sphere  illuminated  by  a  plane  wave. 


2  Fock  Currents  on  a  Convex  Cylinder 

We  now  derive  simple  expressions  for  the  Fock  currents  [12]  on  a  convex  perfectly  conducting 
cylinder  illuminated  by  an  obliquely  incident  plane  wave  (see  Figure  1).  These  expressions 
are  determined  by  applying  the  procedure  in  [13]  to  the  Fock  currents  that  are  valid  for  normal 
incidence.  The  Fock  currents  for  normal  incidence  are  determined  from  expressions  in  [14].  Some 
of  the  results  derived  in  this  section  are  available  in  the  Hterature,  for  example  [15],  [16,  pp.83- 
86],  and  [14].  We  note,  however,  that  the  Fock  currents  for  oblique  incidence  in  the  illuminated 
region  are  given  incorrectly  in  [15]  and  [14],  and  are  not  given  in  [16],  and  there  are  other  errors 
as  well.  In  the  following  derivations  of  the  Fock  current  approximations,  the  cylinder  axis  is 
parallel  to  z,  and  the  unit  normal  vector  for  the  cylinder  is  denoted  n.  denotes  the  radius 
of  curvature  of  the  curve  that  describes  the  cylinder  cross-section  in  a  plane  perpendicular  to  z, 
and  it  is  assumed  that  pn  is  a  continuous  function.  The  incident  plane  wave  propagates  in  the 
direction  k*.  An  time  dependence  is  assumed  throughout  with  a;  >  0. 

2.1  The  Fock  Current  for  TE  Oblique  Incidence 

In  this  case  the  incident  electric  and  magnetic  plane-wave  fields  are  given  by 

E^^-‘(r)  = (1) 
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Figure  1:  Cross  section  of  a  convex  cylinder  with  smoothly  varying  radius  of  curvature. 

and 

=  (2) 

where  Zq  is  the  impedance  of  free  space,  k’  =  and  k’  =  -x  sin  cos  4>^  -y  sin  0*  sin  <j>‘  -z  cos  0" 

^  i  ^  i 

is  the  propagation  direction  of  the  incident  plane  wave.  Moreover,  9  and  <f>  are  the  spherical 
unit  vectors  corresponding  to  the  direction  (0%^®).  We  let  r  =  xx  +  yy  +  zz  denote  a  general 
point  and  let  Tq  =  xx  +  yy  denote  a  point  in  the  plane  z  =  0. 

In  the  special  case  of  normal  incidence  (0®  =^90°)  all  fields  are  independent  of  2;  and  the 
direction  of  propagation  of  the  incident  field  is  k*  =  k®’®®  =  -xcos(^®  -  ysin<^®,  so  that  the 
incident  magnetic  field  is  given  by  H^^’®’®®(ro)  =  ze®^'’"'*^® .  The  Fock  current  in  the  illuminated 
region  can  be  written  a.s  [14,  p.l666] 

=  n  X  H^^’®'®®(ro)  G  ^  n  •  k®-"  j  (3) 

where  G{x)  =  e®^/^£r(x),  with  g{x)  being  the  Fock  function  defined  in  [17,  pp.  63-64].  These  Fock 
functions  (as  well  as  f{x)  and  F{x)  defined  below)  are  most  easily  obtained  by  interpolation  of 
tabulated  values  for  jxj  <  3,  and  by  use  of  asymptotic  expressions  for  |x|  >  3.  Note  that 
n  X  H^^-®’®®(ro)  =  where  is  the  tangent  unit  vector  that  is  perpendicular  to  z  and 

satisfies  t„  =  n  x  z.  Since  n  x  H^^’®’"(ro)  is  simply  one  half  of  the  PO  current  at  the  point  Tq, 
the  Fock  current  (3)  can  also  be  written  as 

=  1  .  k-"j  .  (4) 
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Equation  (4),  valid  for  normal  incidence,  will  now  be  used  in  conjunction  with  the  procedure 
in  [13]  to  obtain  the  Fock  current  for  oblique  incidence.  According  to  [13,  Eq.(13)],  for  oblique 
incidence  the  current  can  be  expressed  as 


K^'^(r)  =  in  sin  6'  +  z 


i  cos 
k  sin 


—ikz  cos  6^ 


J  k’^ksm  6^ 


(5) 


where 


/C”(ro)], 


ifc  fe  1x6'  ^  must  be  replaced  by  A;sin0'  in  the  expression  for 

In  [13]  it  is  also  shown  that  applying  (5)  to  the  PO  current  for  normal  incidence 
yields  the  PO  current  for  oblique  incidence.  Thus,  by  applying  (5)  to  (4)  we  obtain  the  following 
expression  for  the  Fock  current  for  oblique  incidence 


At 


G\ 


kpn{Yo) 


-|l/3 


nk* 


-I  k-^ksmB^ 


+  z 


i  cos  0® 


2k  sin  O'- 


T  e 


—ikz  cos  0^ 


^TE,PO,n^ 

^  ^^dtn 


d 


kpniro) 


-|l/3 


n  •  k®-” 


(6) 


A: —►A;  sin 


where  =  2  h  x  H^^’®(r)  is  the  PO  current  for  oblique  incidence.  Using  the  relation 

•  k®’®®  =  k®’®®  •  in/ Pni  one  finds  that  for  large  kpn  the  second  term  of  (6)  is  negligible  compared 
to  the  first  term  of  (6).  Since  sin0®  n  •  k®’®®  =  n  •  k®,  we  find  that 


'kpn{ro)\ 

.  2  ) 


1/3 


n  •  k®-®® 


'kpg{r)' 


1/3 


nk® 


where 


J  k^k  sin0* 


^  Pni'^o) 

^ 


(7) 


(8) 


is  the  radius  of  curvature  of  the  geodesic  that  forms  the  angle  0®  with  the  z  axis.  Consequently, 
in  the  illuminated  region,  the  Fock  current  for  oblique  incidence  is  given  by 


K^®(r)  =  n  X  H^^-®(r)  G 


kPgjr) 

2 


11/3 


nk® 


(9) 


The  total  current  approximated  in  (9)  is  the  sum  of  the  PO  current  and  the  nonuniform  (NU) 
current.  By  defining  the  modified  Fock  function 

G{x)  =  G{x)  -  2  (10) 

and  using  G  instead  of  G  in  (9)  we  obtain  the  Fock  current  approximation  for  the  TE  NU  current 
in  the  illuminated  region 


K^^-^^(r)  =  n  X  H^^’®(r)  G 


(11) 
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Since  fi  •  k*  =  0  everywhere  on  the  shadow  boundary,  the  Fock  current  at  the  point  Vsb  on  the 
shadow  boundary  is 

K^^irsb)  =  nsb  X  (12) 

The  asymptotic  formula 


G{x)  =  2  +  ^  +  0{x  ®),  X  -oo 

shows  that  in  the  illuminated  region  far  from  the  shadow  boundary  the  Fock  current  is 

1  \ 


TP-/  V  A  t-ttp’,/  ^  •  n  X  H^'^’Tr)  _ 

r  =  2  n  X  ^ ^  +  0 

kpg{r){n  ■ 


(13) 


(14) 


In  (14)  the  first  term  is  the  PO  current  and  the  second  term  is  the  first  higher-order  optics  term. 

Now  consider  the  Fock  current  on  the  shadow  side  of  the  cylinder.  For  normal  incidence  the 
formulas  of  [14,  p.  1666]  show  that  the  Fock  current  is 

K^^'”(ro)  =  Fon,i  -  E^^-*>(r,6o)  (15) 


where  Fsto  is  the  point  of  the  shadow  boundary  that  hes  in  the  z  =  0  plane,  Yq  =  1/Zo  is  the 
admittance  of  free  space,  ns6  is  the  normal  at  the  shadow  boundary,  Sn(ro)  is  the  distance  along 
the  cylinder  from  the  point  r^bo  on  the  shadow  boundary  to  the  point  Fq,  and 


""■^27  /0n(s„)2/3’ 


(16) 


A  ^ 

is  the  Fock  parameter  (the  integration  in  (16)  is  on  the  surface  of  the  cylinder).  Moreover,  t„  is 
that  unit  tangent  vector  to  the  cylinder  which  is  perpendicular  to  z  and  is  equal  to^  k*’”  at  the 
shadow  boundary.  On  the  illuminated  side  the  unit  tangent  vector  was  defined  as  tn  =  h  x  z, 
and  the  two  unit  tangent  vectors  on  either  side  of  the  shadow  boundary  are  therefore  not  always 
equal  at  the  shadow  boundary  (in  some  cases  they  are  opposite).  However,  the  unit  tangent 
vector  in  the  illuminated  region  was  used  only  in  intermediate  steps  of  the  derivation,  and  it 
does  not  appear  in  any  of  the  final  formulas  for  the  current  in  the  illuminated  region.  Note  also 
that  Yq  Usb  •  E^^’*’”(rs6o)  =  since  the  incident  electric  field  is  always  parallel  to  the 

normal  at  the  shadow  boundary. 

To  obtain  the  Fock  current  for  oblique  incidence,  we  insert  (15)  into  (5)  and  note  that  (5) 
remains  unchanged  when  t„  is  reversed.  It  is  found  that  for  large  kpn,  the  terms  involving  the 
derivative  of  g{xn)  or  the  derivative  of  [/3n(rs6o)//3n(ro)]^'^^  are  negligible  compared  to  the  two 
other  terms  that  result  when  (15)  is  inserted  into  (5).  Thus,  we  have 


K^£^(r)  =  Fo 


cin  Jksine'sniTo)  ,  a£cOS^  6  iksmB'sn(To) 

t^sine^e  +^ksm6^  dtj 


(17) 
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The  last  factor  in  (17)  simply  equals  (t„sin^*  —  where 

is  the  tangent  along  the  geodesic  that  forms  the  angle  with  the  2:  axis,  as  shown  in  Figure  1. 
By  unrolling  the  cylinder  one  finds  that  sin^’5„(ro)  -  cos  9\z  -  -  sair),  where  S5(r)  is  the 

distance  along  the  geodesic  from  the  point  =  Ss{,x  +  j/^jy  +  on  the  shadow  boundary  to 
the  point  r  in  the  shadow  region.  Thus, 

g-ii:2cose'  Ltfcs„(ro)  gik''"-rs6ol  _  g»fc«s(r)  gik'-r^b  /■.o') 

1  Jfc-*fcsine‘ 

showing  that 


„—{kz  COS  9 


‘  [n,, .  ^  =  na  ■ 

j  aC""^ /c  sm  w 


iksg(r) 


(19) 


From  the  relation  (8)  between  and  pg  and  the  formula  =  sin^Vsg  we  find  that  [x„]fc_* 
simply  equals 


.  =  ('iV'’  r 

\2y  </r,6  Pgisgf!^ 


ksinO^ 


(20) 


where  the  integration  is  performed  along  the  geodesic  that  starts  at  rgb  and  ends  at  the  point 
r  in  the  shadow  region.  Inserting  all  these  results  into  (17)  proves  that  the  Fock  current  in  the 
shadow  region  is  given  by 


K^^(r)  =  Fo  nsb  ■  t,  (21) 

when  the  incident  field  is  the  plane  wave  in  (1)  and  (2). 

At  the  shadow  boundary  the  expression  in  (21)  reduces  to 

K^^(r,6)  =  Yoh,b  ■  E^®’‘(r,,)  g(0)  k*  (22) 

which  can  be  shown  to  be  in  agreement  with  (12)  by  use  of  ^(0)  =  G(0),  =  Asb  or  =  —6*6, 

and  tg^sb  =  k‘.  For  large  positive  values  of  its  argument  the  Fock  function  g{x)  is  given  by  [17, 
p.  64] 

sW~ftS(_ft)-  (23) 

where  is  the  negative  of  the  first  zero  of  the  derivative  of  the  Airy  function.  Substituting  (23) 
into  (21)  gives  us  the  first  creeping  wave  on  the  convex  cylinder 


K^^(r)  ~  Fo  n*6  •  E^^-'(r*6) 


/?iAz‘(-^i) 


(  P9i^sb)\ 

\  Pgir)  ) 


1/6 


''9- 


(24) 


6 


2.2  The  Fock  Current  for  TM  Oblique  Incidence 

In  this  case  the  incident  electric  and  naagnetic  plane- wave  fields  are  given  by 


and 


E™’’(r)  =  -e' 
H™-’(r)  =  Fo 


(25) 

(26) 


where  Yq  is  the  impedance  of  free  space.  For  normal  incidence  the  electric  field  in  (25)  is  given  by 
E^^’®’"(ro)  =  and  [14,  p.l666]  shows  that  the  Fock  approximation  to  the  total  current 

on  the  illuminated  side  of  the  cylinder  is 


i!:J^’”(ro)  =  iYoi-  E^^’‘’"(ro)  F 


kpuif^o) 


1/3 


n  •  k*’” 


kpn{To) 


-1/3 


(27) 


where  F{x)  —  e*®*/^/(a:)  and  f[x)  is  the  Fock  function  defined  in  [17,  pp.  63-64].  Since  the  PO 
current  for  normal  incidence  is 


K™'PO'-^{vo)  =  2  fi  X  H™’*'"(ro)  =  -2Fon  •  k*’”  z 

we  can  write  the  Fock  current  (27)  as 
K^^’"(ro)  =  zK™'^{ro) 

'kpn{roy^^^ 


(28) 


kpniro) 


■11/3 


nk 


.i,n 


-1 


nk* 


.(29) 


According  to  [13]  the  current  for  oblique  incidence  can  be  determined  from  the  current  for 
normal  incidence  through  the  equation 


K™{r)  =  e 


—ikz  COS  6' 


K™-’{to)] 


k^ksinO* 


(30) 


and  in  [13]  it  is  also  shown  that  applying  (30)  to  the  PO  current  for  normal  incidence  yields  the 
PO  current  for  oblique  incidence.  Inserting  the  Fock  current  (29)  into  (30)  shows  that  the  Fock 
current  for  oblique  incidence  is 

K^^(r)  =  zFrJ^(r)  =  -inx  H^^-*(r) 

where  we  have  used  (7)  and  inserted  the  expression  K™’^‘^(r)  =  2  n  x  H^^’’(r)  for  the  PO 
current.  Alternatively,  we  can  insert  the  Fock  current  (27)  into  (30)  to  get 


1/3 


n-k* 


kp^ 


■11/3 


-1 


nk* 


(31) 


iYo 


sin^  9^ 


z  ■  E™-*(r)  F 


kPgjr) 

2 


2 


-1/3 


(32) 
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where  we  have  used  the  formula 


^—ikz  cos 


[z  • 


k—^ksiuB 


z  •  E^^’»(r) 
sin  0^ 


Since  the  total  current  approximated  in  (31)  or  (32)  is  the  sum  of  the  PO  current  and  the 
NU  current,  by  defining  the  modified  Fock  function 

F{x)  =  F{x)  —  2ix  (34) 

and  using  F  instead  of  F  in  (31)  or  (32)  we  obtain  the  Fock  current  approximation  for  the  TM 
NU  current  in  the  illuminated  region 


K™-^^(r)  =  z/i:™-^^(r) 


=  n  X  H^^-’(r)  n  •  k* j  f 


2 


KTM,nu^^ 


sin^  9^ 


z  •  E™’*(r)  F  f  '  n  •  k* 


2 


At  a  point  r^b  on  the  shadow  boundary  both  n  •  k’  and  n  x  equal  zero,  and  we  find 


UJ-LU/U 

n  X  H^^’»(r)  yp  ^  _  h^b  •  H^^’»(r,fe) 

n-k*  sin0‘  n^t-^'sin^* 

so  that  the  Fock  current  (31)  at  the  shadow  boundary  becomes 


z,  r  Tsb 


K™{rsb)  =  zK™{rsb)  =  z 


F(0) 

Usb  •  ^  sin  6^  \  ^ 


which  is  seen  from  the  plane-wave  expressions  (25)  and  (26)  to  agree  with  the  result  obtained 
from  (32). 

For  laxge  negative  values  of  its  argument,  the  function  F{x)  is  given  by  [17,  p.  64] 


F{x)  ~  2ix  +  —  0{x  ®), 


x  — oo. 


so  that  the  current  in  the  illuminated  region  away  from  the  shadow  boundary  is 
K^CD  ~  2  n  X  +  0  1  ^ 

kpg{r)  •  k‘)  V(^Pir)^(n  •  k‘)®> 


in  which  the  first  term  is  the  PO  current  and  the  second  term  is  the  first  higher-order  optics 
term. 
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On  the  shadow  side  of  the  cylinder,  the  Fock  approximation  for  normal  incidence  to  the  total 
current  is  [14,  p.  1666] 

A-™."(ro)  =  i  /(x„)e“-(»>  (41) 

Asb-^  V  Pn(^o)  J  \  ^  J 

where  the  Fock  parameter  is  defined  in  (16).  Inserting  (41)  into  (30),  and  using  (18),  we  find 
that  the  Fock  current  for  oblique  incidence  is 


(r)  Ash  •  H^^’^'(r,6)  / Pg(r^6)y^^  f  kpgjr) 

Asb  ■  sinO^  \  Pai^)  /  V  2 


-1/3 


(42) 


with  Xg  given  by  (20). 

At  the  shadow  boundary,  the  value  of  the  Fock  current  (42)  is 


A:™(r,,)  =  i/(0) 


ns6  •  H™’’(rs6)  ( kpgijCsb)' 


-1/3 


hsb  •  sin  (9* 


(43) 


which  agrees  with  the  result  in  (38)  because  /(O)  =  F'(O).  For  large  positive  values  of  its 
argument,  the  Fock  function  f{x)  is  given  by  [17,  p.  64] 


/(^) 


X  — )■  +OC, 


(44) 


where  ai  is  the  negative  of  the  first  zero  of  the  Airy  function.  Substituting  (44)  into  (42)  gives 
the  first  creeping  wave  on  the  convex  cylinder 

i  n.,  ■  f  .  (45) 


Ai'(— q;i) 


Asb  ■  ^  sin  V  Pgiy) 


2.3  Examples 

We  now  specialize  the  Fock  current  approximations  for  a  general  convex  PEC  cylinder  obtained 
above  to  the  cases  of  a  circular  and  parabolic  cylinder.  Consider  first  a  circular  cylinder  of  radius 
a  illuminated  by  the  TE  plane  wave  given  by  (1)  and  (2).  For  the  Fock  current  approximation 
to  the  TE  NU  current  in  the  illuminated  region  we  obtain  from  (11) 


K^'®’^^(a,  4>,  z)  ~  sin  +  z  cos  9^  sin(^  —  <j)'') 


ifca  sin  0*  cos((?{>— —ikz  cos  9^ 


■G 


'  ka  sin  9 


.•\  1/3 


cos(<^  —  <f>^) 


<n/2  (46) 


9 


1.6 


Figure  2:  Amplitude  of  {a,(i),z)  for  circular  cylinder,  ka  =  60,  =  180°,  9^  =  90°. 


while  for  the  shadow  region  (21)  yields 

K^^(a,  z)  Ri  [-^  sin  ±  z  cos  0*]  g-»fc2cose>' 


^fcasin  9^ 


1/3 


|(?i^'±7r/2-<^| 


^/2  ^  1^  —  ^  TT 


(47) 


where  the  i  corresponds  to  the  shadow  boundary  at  <^‘  respectively.  For  illumination  by 

the  TM  plane  wave  given  by  (25)  and  (26),  for  the  illuminated  region  we  obtain  from  (35) 


K™’'^^(g,  (?i,  .^)  R!  izFo  e-ikzcosB' 


F  —  ^Arasin^’j  ^  cos(<^— ^®) 


(|A;asin^®)  ^ 

|^-^*|<7r/2  (48) 


while  for  the  shadow  region  (42)  yields 

Kj^{a,(j),z)  Ri  iYo  ^-ikzcose 


.  /  ^|A;asin^’j  ^  |(/>‘ ±  7r/2  —  (^| 

^|A:osin^*^  ^ 

7r/2<  <x. 


(49) 


In  Figure  2  we  have  plotted  the  amplitude  of  the  approximate  TE  NU  current, 
j^TE,NU^a^  ^  circular  cylinder  with  ka  =  60  obtained  from  (46)  and  (47)  with  (f>^  =  180° 

and  9"  =  90°  (normal  incidence)  along  with  the  amplitude  of  the  exact  current  obtained  from  the 
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<f>  (degrees) 

Figure  3:  Amplitude  of  for  a  circular  cylinder,  ka  =  60,  =  180°,  =  30°. 

eigenfunction  expansion  for  scattering  from  a  circular  cylinder.  The  amplitude  is  discontinuous 
at  the  shadow  boundary  a.t  (f>  =  90°  because  the  TE  PO  current  is  discontinuous  there,  and 
hence  when  subtracted  from  the  continuous  total  current  gives  a  discontinuous  NU  current.  The 
approximate  current  is  seen  to  closely  fit  the  exact  current  except  for  the  ripple  behavior  of  the 
exact  current  in  the  deep  shadow  region  in  the  vicinity  of  (^  =  0°.  This  ripple  behavior  is  the 
result  of  interference  between  the  dominant  creeping  waves  launched  aX  <f)  =  ±90°  both  of  which 
are  present  in  the  exact  eigenfunction  expansion  of  the  current,  whereas  the  approximate  current 
in  the  range  0°  <  <f>  <  180°  contains  only  the  dominant  creeping  wave  launched  at  <f>  =  90°. 
Much  the  same  remarks  apply  to  the  comparison  of  the  plots  in  Figure  3  of  the  approximate  and 
exact  z  component  of  the  TE  NU  current  for  oblique  incidence  (^  =  30  ). 

In  Figure  4  we  have  plotted  the  amplitude  of  the  approximate  TM  NU  current, 
obtained  from  (48)  and  (49)  for  normal  incidence,  along  with  the  amplitude 
of  the  exact  current.  Here  there  is  no  discontinuity  at  the  shadow  boundary  because  the  TM 
PO  current  is  zero  in  the  illuminated  region  at  the  shadow  boundary.  The  approximate  cur¬ 
rent  closely  matches  the  exact  current  everywhere.  No  ripple  behavior  is  observed  here  for  the 
exact  current  because  the  creeping  waves  attenuate  much  more  rapidly  for  TM  than  for  TE 
illumination. 

Next  we  consider  the  parabolic  cylinder  with  focal  length  h.  The  currents  are  given  in 
terms  of  the  parabolic  cylindrical  coordinates  (^,»?,^)  discussed  in  [17,  ch.  7].  The  surface  of 
the  scattering  cylinder  is  given  hy  rj  =  rji  =  y/^.  The  derivations  of  the  Fock  currents  are 
straightforward  but  somewhat  lengthy.  Expressions  given  in  [17,  ch.  7]  are  used,  as  well  as  the 
formula  pn  =  (^^  ±  for  the  radius  of  curvature  of  a  parabolic  cylinder  in  the  plane 

normal  to  the  z-axis.  We  give  only  the  results  here.  For  illumination  by  the  TE  plane  wave  given 
by  (1)  and  (2)  with  0  <  (f)^  <  tt,  we  obtain  from  (11)  the  Fock  approximation  for  the  NU  current 
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(f>  (degrees) 

Figure  4:  Amplitude  of  (a,  4>i^)  for  a  circular  cylinder,  ka  —  60,  =  180°,  =  90°. 


in  the  illuminated  region 


7/1,2:)  «  l^sin^*  +  zcos^ 


.  f  2  __2 

COS  (j)^  +  Tji  sin  (f>‘^  —ikh  sin  6'  (cos0*- — 5-^+2  sin 


ie + \ 

J,  •  oA  1/3  7/1  cos  sin 


,^-ikz cosei  g  (Msin^y  ^  ^  ,  ^>7/icotf.  (50) 

For  the  shadow  region,  (21)  yields 
[«sine-  +  zcose-] 

■g  fasiny)*^°li'(’?r^(-{  + V'{^  +  ’?iO)t'^’)]  yl  (-;^  i)  .  {  <  Iji cot cA*.  (■ 


•p  sin  In  [t]^  {-(  +  +  r](  cot  J  ^ 

It  can  be  convenient  to  replace  the  variable  ^  which  ranges  from  — oo  to  +oo  with  a  generalized 
angle  variable  (not  to  be  confused  with  the  polar  angle  (f>)  defined  by  (f>g  —  (or  ^  = 

i]i  tan  (j)g)  which  ranges  firom  — 7r/2to7r/2  with  the  shadow  boundary  given  by  (f>g,sB  =  7r/2  — 

In  terms  of  4>g  (50)  and  (51)  become 


kh,  z)  PS  sin  0*  +  z  cos  cos  ^p|(sin  (j)'  +  cos  4>'  tan  <j)g^ 
.g-ifcftsme‘(cos^‘(tan2^-l)+2sm^'tan^,)  ^-ikzcosr  Q  sin  0*)^ (cOS  (f>^  -  sin  ?iHan  , 


<f>g>Tl2-<j>^  (52) 


K^^((^5,A:/i,2)  fs  flsin^' -f  zcos^'l  i)°°‘(<^V2)]-taii^p|cos0g|  1} 
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—ikz  cos  9' 

•e  y 


[kh  sin^’^  ^  In  tan  (f)g  +  |  cos  <j)g\  cot 


cos 


-v/sin0®  ’ 
(t>g<w/2-cl)\  (53) 


For  illumination  by  the  TM  plane  wave  given  by  (25)  and  (26),  in  the  illuminated  region  we 
obtain  from  (35)  or  (36) 


Ri  iYo- 


lin  B' 


—ikh  sin  (  cos 


+2  sin<^ 


_—ikz  cos  i 


ie  +  rilY'^ 

sB'  ^  {rji  cos  (f)^  —  ^  sin  ^*)] 

(kh  cos 

while  in  the  shadow  region  we  obtain  from  (42) 

7/1,2)  «  iYo 


’71/ 


,  ^  >  771  cot 


(54) 


^—ikz  cos  0* 


/  (kh  sin  (9’)^^^  In  (rj^  ^(-(  +  y/^^  +  ?7i  cot  f  ^ 
\/sin  <f>Ykh  sin  0*)^/^ 


2  \  3/4 

,  ^<»7iCot<;ih  (55) 


In  terms  of  the  generahzed  angle  variable  ^g,  these  expressions  are  respectively 

/i:J^-^^((^^,;bh,z)  «  ^-ikzcosB' 

F  ^{kh  sin  9^^^^  (cos  —  sin  tan  <j>g)^ 

(khsin.  0*)V3 


-,  4t,  >  7r/2  -  f 


(56) 


and 


gifc/i sine* {lii[(-taii.^p+|cos.^j,|  i) cot(.^72)]-tan(^5|cosi^3| 

gt  f  [(fc/isin^®)^'^^ln  [(—  tan  <l>g  +  |  cos  4>g\~^)  cot  ^®]j 


^ikz cos  9^ 


■\/sin  <f>Ykh  sin  0’)^/^ 


,  c^g  <  7r/2  -  (57) 


In  Figure  5  we  show  plots  of  the  approximate  TE  NU  current  {4>g,  kh,  z)  for  a  parabolic 

cylinder  with  kh  =  60,  obtained  from  (52)  and  (53)  with  (Y  =  30°  and  0®  =  90°  (normal 
incidence).  Also  plotted  is  the  exact  current  obtained  from  [17,  Eq.(7.47)]  using  expansions 
given  in  [18]  and  [19].  The  approximation,  while  not  quite  as  good  as  for  the  circular  cylinder, 
is  still  very  close. 

In  Figure  6  we  show  a  comparison  of  the  z  component  of  the  approximate  TE  NU  current 
for  the  same  parabolic  cylinder  with  =  30°  and  9^  =  30°  (oblique  incidence),  and  the  exact 
current  obtained  from  [17,  Eq.(7.47)]  together  with  [13],  [18],  and  [19].  Again  the  approximation 
is  quite  good  throughout  the  range. 

Plots  of  the  approximate  TM  NU  current  {(f>g,kh,  z)  obtained  from  (56)  and  (57)  for 

normal  incidence  and  of  the  corresponding  exact  current  obtained  from  [17,  Eq.(7.12)]  and  [18] 
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Figure  7:  Amplitude  of  kh,z)  for  a  parabolic  cylinder,  kh  =  60,  <j>^  =  30°,  0^  -  90°. 

are  displayed  in  Figure  7.  Here  the  approximation,  while  acceptable,  differs  somewhat  from  the 
exact  current  especially  in  the  shadow  region,  primarily  because  the  value  of  the  TM  current 
at  the  shadow  boundary  4>9  =  60°  is  significantly  less  than  the  exact  value.  The  reason  for 
this  discrepancy  is  that  the  Fock  current  approximation  is  derived  under  the  assumption  that 
the  derivative  of  the  radius  of  curvature  at  the  shadow  boundary  is  equal  to  zero,  a  condition 
satisfied  by  the  circular  cylinder  but  not  by  the  parabolic  cyhnder  except  when  the  shadow 
boundary  is  at  the  vertex  of  the  parabolic  cylinder.  The  discrepancy  between  the  exact  and 
approximate  values  at  the  shadow  boundary  of  the  parabolic  cylinder  is  present  in  Figures  5  and 
6  as  well,  but  to  a  lesser  degree. 


3  Integrating  the  Fock  Currents  to  Obtain  ILDC’s 

We  will  now  show  how  the  ILDC’s  corresponding  to  the  Fock  currents  presented  in  Section  2 
can  be  obtained  efficiently.  The  ILDC’s  for  the  shadow  boundary  of  a  convex  cylinder  represent 
the  fields  radiated  by  the  nonuniform  Fock  currents  on  that  cylinder.  Specifically,  the  far  field 
radiated  by  a  strip  of  thickness  dz',  which  is  described  by  the  curve  C  on  the  cylinder,  is  given 

by 

dE^^(r)  y  [K^^{0\(j)\s)-Tr-K^^{e\<t>\s)]  ds,  (58) 

47rr  6  ^ 

where  r  is  the  direction  to  the  far-field  observation  point,  tc{s)  is  the  parametric  representation 
of  C,  and  is  the  nonuniform  Fock  current.  In  the  illuminated  region,  the  curve  C  is  normal 

to  the  shadow  boundary.  In  the  shadow  region,  C  is  the  geodesic  whose  angle  with  the  axia 
direction  is  6\  as  shown  in  Figure  1. 
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In  general  it  is  not  possible  to  evaluate  the  integral  (58)  in  closed  form.  However,  as  we 
will  now  demonstrate,  an  accurate  and  rapidly  calculable  approximation  to  this  integral  can  be 
obtained. 

3.1  Quadratic  Polynomial  Approximations  of  the  Amplitude  and 
Phase  of  the  ILDC  Radiation  Integral 

Begin  by  considering  one  of  the  components  of  the  electric  field  dE  in  (58)  and  write  this  com¬ 
ponent  as 

dE(r)=  r  A{s)e-^^^^Us,  (59) 

5a 

where  the  amplitude  ^(s)  and  the  phase  P(s)  are  functions  of  the  angles  of  incidence  and 
observation.  The  amplitude  and  phase  functions  are  determined  from  the  nonuniform  current 
and  the  far- field  Green’s  function  occurring  in  (58).  For  convenience,  we  allow  the  amplitude 
A(s)  to  be  complex.  Moreover,  the  expressions  in  Section  2  for  the  Fock  currents  and  the 
factor  determine  the  unwrapped  phase  of  the  integrand,  so  that  the  function  P(s) 

is  continuous  except  possibly  at  the  shadow  boundary.  Thus,  the  phase  function  P(s)  is  not 
restricted  to  a  finite  interval  of  length  27r. 

We  can  divide  the  integration  in  (59)  into  subintervals  and  have 


where  Si  =  -Sn+i  =  and  Sj  <  Sj^i.  The  form  of  the  functions  .4(s)  and  P{s)  leads  one  to 
assume  that  these  functions  can  be  accurately  approximated  by  quadratic  polynomials  of  s  in 
each  of  the  intervals  Sj  <  s  <.  when  these  intervals  are  chosen  appropriately.  Hence, 

I  ~  r*'  (cos^  +  cxs  +  (61) 

J Sj  J Sj 

where  Cq,  Ci,  C2,  C3,  C4,  and  C5  are  constants  independent  of  s. 

These  constants  will  be  determined  such  that  the  quadratic  polynomials  fit  A  and  P  at  the 
endpoints  Sj  and  Sj+i  and  at  the  midpoint  5j+i/2  =  (sj  +  Sj+i)/2: 

A(s)  =  CqS^  -t-  CiS  +  C2,  for  S  =  Sj,  S  =  Sj+1/2,  5  =  Sj+i, 

which  gives 

C-2[S,»+S,)  +*  (3,^.  -  3,03 

_  "4(5j+i)  -  A(sj) 

Cl  —  -  -  Co(5j+i  +  5j), 


(62) 

(63) 

(64) 
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Figure  8:  Circular  cylinder  illuminated  by  a  TE  plane  wave 


and 


C2 


- r  CoSjSj+i 


(65) 


Similarly, 

P{s)  =  C3S^  +  C4S  +  C5,  for  5  =  Sj,  s  =  Sj+1/2,  -s  =  Sj+i  (66) 

and  we  find  that  C3,  C4,  and  C5  are  given  by  (63),  (64),  and  (65),  respectively,  with  A  replaced  by 
P.  These  expressions  for  the  constants  co,  Ci,  C2,  C3,  C4,  and  C5  are  well  behaved  for  all  functions 
A  and  P  that  are  continuous  in  the  interval  sj  <  s  <  Sj+i.  In  principle  we  could  let  P(s)  be 
continuous  for  all  s  except  at  the  shadow  boundary.  However,  in  practical  calculations  we  permit 
P{s)  to  have  discontinuities  that  are  multiples  of  27r  at  the  points  Sj. 

In  the  appendix  we  have  shown  how  the  integral  on  the  right  side  of  (61)  can  be  evaluated  in 
closed  form  in  terms  of  the  Fresnel  function.  Thus  the  ILDC’s  are  expressed  as  a  finite  summation 
involving  the  Fresnel  function  and  the  values  of  the  amplitude  and  the  unwrapped  phase  of  the 
integrand  in  (58)  evaluated  at  the  points  Sj. 


3.2  Numerical  Verification  of  the  Quadratic  Approximations 

We  will  now  use  a  numerical  example  to  verify  the  quadratic  approximation  of  the  phase  and 
amplitude  in  (61).  We  consider  the  far  field  radiated  by  the  NU  current  on  the  circular  cylinder 
shown  in  Figure  8.  The  incident  field  is  a  plane  TE  wave  with  HI  =  and  only  the  NU 
current  in  the  angular  region  0  ^  ^  ^  tt  will  be  taken  into  account.  The  magnetic  far  field  for 
this  NU  current  can  be  written  as  Pf  “^’^^(r)  =  j Vkr ,  where  the  far-held  pattern 

is  given  by 

^  n  A{s)e-^^^^Us,  (67) 

Jsa 
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1.6 


0  (degrees) 

Figure  9:  Amplitude  of  the  exact  nonuniform  TE  current  on  circular  cylinder  for  ka  =  60. 
with  Sa  =  0,  56  =  TT, 

^  cos(<^  -  s),  (68) 

and 

P{s)  =  ka  cos(^  —  s)  —  Phase(/ir^’^’^^(s))  (69) 

where  Phase(ii['^^’^^(s))  is  continuous  (i.e.,  not  modulo  27r)  apart  from  the  discontinuity  at  the 
shadow  boundary. 

The  amplitude  and  unwrapped  phase  of  NU  current  are  shown  in  Figure  9  and  10,  respectively, 
for  A;a  =  60.  Notice  that  both  the  unwrapped  phase  and  the  amphtude  of  the  NU  current  have 
discontinuities  at  the  shadow  boundary,  and  that  the  unwrapped  phase  is  almost  linear  on  the 
shadow  side. 

Figures  11-14  show  the  exact  values  and  the  quadratic  approximations  of  A(s)  and  P{s)  for 
4>  =  Q°  and  (j>  =  90°.  The  number  of  segments  n  in  (60)  is  6  and  Sj  =  —  l)/6,  so  the 

shadow  region  0  <  s  <  7r/2,  as  well  as  the  illuminated  region  7r/2  <  s  <  tt,  are  divided  into 
three  subregions.  Thus,  the  length  of  the  subintervals  is  5A.  To  better  illustrate  the  accuracy 
of  the  quadratic  phase  approximation,  we  have  introduced  discontinuities  in  the  phase  at  the 
endpoints  Sj.  We  see  that  the  quadratic  phase  approximation  is  excellent  for  both  these  angles 
of  observation,  but  that  the  quadratic  amplitude  approximation  deviates  somewhat  from  the 
exact  amplitude  for  =  0°. 

The  quadratic  approximations  can  be  improved  by  using  unequal  spacing  with  a  greater 
density  of  points  taken  in  the  region  of  s  that  contributes  most  strongly  to  the  far  field  at  <f>. 
As  an  example,  this  is  done  in  Figures  15  and  16  for  ^  =  0°  where  the  endpoints  are  given  by 
si  =  0,  52  =  60°,  53  =  80°,  54  =  90°,  55  =  105°,  53  =  135°,  and  57  =  180°.  The  quadratic 
amplitude  approximation  is  clearly  better  in  Figure  15  than  in  Figure  11  with  equal  spacing. 
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EXACT 

n=6 


Figure  12:  Exact  values  and  quadratic 
The  endpoints  are  sj  =  7r{j  —  l)/6. 


approximation  of  P(s)  (or  n  =  6,  <f>  =  0^  and  Jru  =  60 


Figure  13:  Exact  values  and  quadratic  approximation  of  /Ifa)  for 
The  endpoints  are  sj  =  ir(j  —  l)/6. 


n  —  Q,  <f>  =  90®,  and  ka  =  60. 


2l 


J] 


Figure  14:  Exact  values  and  quadratic  approximation  of  P(s)  for  n  =  6,  ^  =  90*’,  and  ka  =  60. 
The  endpoints  are  Sj  =  7r(j  —  l)/6. 


Figure  15:  Exact  values  and  quadratic  approximation  of  A(s)  for  n  =  6,  ^  =  0°,  and  ka  —  60. 
The  endpoints  are  Si  =  0,  52  =  00°,  53  =  80°,  S4  =  90°,  5$  =  105°,  s&  —  135°,  and  57  =  180°. 


Figure  16:  Exact  values  and  quadratic  approximation  of  P{s)  for  n  =  6,  —  0°,  and  ka 

The  endpoints  are  =  0  ,  52  =  60°,  53  =  80°,  54  =  90°,  55  =  105°,  se  —  135  ,  and  57  —  181 


0  100  200  300 


(f)  (degrees) 

Figure  17:  Exact  and  approximate  NU  TE  far-held  pattern  for  ka  =  60  and  n  =  4,  n  = 
n  =  8.  The  endpoints  are  Sj  =  Tr(j  —  l)fn]  0°  <  ^  <  360°. 


^  (degrees) 


Figure  18:  Exact  and  approximate  NU  TE  to^Held  pattern  for  ka  =  60  and  n  =  4,  n  -  6.  and 
n  =  8.  The  endpoints  are  Sj  =  7r(i  -  l)/n,  0  <  ^  _  1 

Future  work  could  deterndne  the  opttaal  .t  of  -^poinfo  \aa  a  fonct^  l"prtbl!y 

3  -rr •  i -sir. 

and  p  I,  a  TP  nlanp  wave  incident  on  a  circular  cylinder,  the  JNU 

17  and  IS.  In  th»e  6^res._the  ex.t  NW^f-field  p^n  is 

and  n  =  8,  and  the  endpomts  axe  grven  by  aj  -  xO  )/  _  ^  impossible 

n  =  8,  and  we  conclude  that  for 

^  Thi  iba  r  .h*  mwai...  toXniqu. 

order  of  A/ 6. 
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Figure  19:  Three  dimensional  scatterer  illuminated  by  a  plane  wave. 

4  The  Use  of  the  ILDC’s  for  a  3D  Scatterer 

In  this  section  we  show  how  the  ILDC’s  for  a  convex  cylinder  are  applied  to  calculate  the  scattered 
far  field  of  a  3D  scatterer.  In  applications  the  PO  field  will  generally  be  known,  and  what  is 
desired  is  the  scattered  field  due  to  the  NU  currents  on  the  3D  scatterer  in  the  vicinity  of  its 
shadow  boundary.  The  total  scattered  field  is  then  obtained  by  adding  the  NU  scattered  field  to 
the  PO  field.  Consider  a  3D  scatterer  with  a  smoothly  varying  radius  of  curvature  illuminated 
by  a  plane  wave  (see  Figure  19).  The  illuminated  and  shadowed  portions  of  the  surface  are 
separated  by  the  shadow  boundary  curve  on  the  surface  of  the  scatterer.  Consider  a  point  P 
on  the  shadow  boundary  and  generate  a  curve  C  on  the  surface  of  the  scatterer  through  P 
motivated  by  the  behavior  of  the  currents  on  the  surface.  In  the  shadow  region,  C  is  the  geodesic 
emanating  from  P  in  the  direction  of  the  illuminating  plane  wave  —  the  curve  followed  by  the 
currents  -  while  in  the  illuminated  region  where  the  currents  are  determined  locally,  C  is  simply 
the  intersection  of  the  surface  with  the  plane  through  P  perpendicular  to  the  shadow  boundary. 
A  2D  convex  cylindrical  (in  general  not  circular-cylindrical)  surface  can  then  be  generated  by 
passing  straight  fines  (generators)  through  each  point  on  C  parallel  to  the  shadow  boundary  at 
P.  Now  generate  a  curve  C  through  P  on  the  cylindrical  surface  in  the  same  way  that  C  was 
generated  on  the  surface  of  the  3D  scatterer.  Assume  that  the  radius  of  curvature  of  the  3D 
scatterer  is  large  in  the  direction  parallel  to  the  shadow  boundary  at  P  for  points  on  the  geodesic 
part  of  C,  especially  in  the  vicinity  of  P.  Then  the  geodesic  portion  of  C  in  the  vicinity  of  P 
will  be  closely  approximated  by  the  geodesic  portion  of  C  (a  generalized  helical  curve  on  the 
convex  cylindrical  surface  making  a  constant  angle  with  its  generators),  and  the  ILDC’s  of  the 
cylindrical  surface  at  P  will  be  a  good  approximation  to  the  fields  radiated  by  the  NU  currents 
on  a  differential  strip  of  the  3D  scatterer’s  surface  containing  C.  The  procedure  for  obtaining  the 
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ILDC’s  corresponding  to  the  point  P  on  the  shadow  boundary  of  the  3D  scatterer  is  repeated  for 
all  points  on  the  shadow  boundary,  and  the  ILDC’s  then  integrated  along  the  shadow  boundary 
to  obtain  an  approximation  to  the  far  field  radiated  by  the  NU  currents  on  the  3D  scatterer. 

It  should  be  noted  that  a  creeping  wave  geodesic  on  the  surface  of  the  actual  scatterer  will 
in  general  diverge  from  the  corresponding  hehcal  creeping  wave  geodesic  on  the  approximating 
convex  cylinder  surface  as  the  distance  from  the  launching  shadow  boundary  point  increases. 
However  the  exponential  decay  of  the  creeping  waves  (see  Eqs.  (24)  and  (45))  insures  that  the 
major  contribution  to  the  diffracted  far  field  comes  from  the  current  on  the  portion  of  the  geodesic 
close  to  the  shadow  boundary.  Here  the  generalized  helical  geodesic  can  be  expected  to  be  a 
good  approximation  to  the  actual  creeping  wave  geodesic,  especially,  if  as  assumed,  the  radius  of 
curvature  of  the  actual  3D  scatterer  is  large  in  the  direction  parallel  to  the  shadow  boundary  for 
points  on  the  geodesic  close  to  the  launching  shadow  boundary  point.  An  important  consequence 
of  this  concentration  of  the  creeping  wave  current  in  the  vicinity  of  the  shadow  boundary  is  that 
the  major  ILDC  contribution  to  the  far  field  is  in  the  forward  direction,  since  the  field  in  this 
region  can  be  viewed  as  due  to  the  shedding  of  the  creeping  waves  along  rays  tangential  to  the 
scatterer  at  points  close  to  the  shadow  boundary. 

4.1  Calculation  of  the  Fields  Scattered  by  a  Sphere 

As  an  application  of  the  NU  ILDC’s  for  the  shadow  boundary  of  the  PEC  circular  cylinder,  we 
calculated  the  far  field  of  a  PEC  sphere  illuminated  by  a  plane  wave  by  adding  the  integral  of 
the  ILDC  for  a  cylinder  to  the  PO  field  of  the  sphere,  and  compared  the  result  to  the  exact 
scattered  far  field  of  the  sphere  as  obtained  from  the  Mie  series  solution.  The  sphere,  of  radius 
a,  is  assumed  centered  at  the  origin  of  a  Cartesian  coordinate  system  and  illuminated  by  a  plane 
wave  propagating  from  above  in  the  negative  z-direction.  The  shadow  boundary  of  the  sphere  is 
thus  a  circle  of  radius  a  in  the  z  =  0  plane.  A  circular  cylinder  of  radius  a,  enclosing  the  sphere, 
whose  axis  lies  in  the  plane  z  =  0,  will  be  tangent  to  the  sphere  along  the  circle  of  radius  a, 
normal  to  the  cylinder  axis,  with  its  center  at  the  origin.  Hence,  a  strip  of  the  sphere’s  surface 
of  constant  equal  width  on  either  side  of  the  meridian  defined  by  ^  =  constant  coincides  in  the 
limit,  as  the  width  of  the  strip  shrinks  to  zero,  with  an  azimuthal  differential  strip  of  surface 
of  the  enclosing  cylinder.  As  the  cyhnder  is  rotated  around  the  z-axis  the  entire  surface  of  the 
sphere  is  obtained.  Thus  the  far  field  of  the  sphere  can  be  obtained  by  adding  the  integral  of  the 
NU-current  ILDC  of  the  cylinder  to  the  PO  far  field  of  the  sphere,  as  the  cylinder  axis  is  rotated 
from  (f>  =  0  to  (j)  =  27r.  (It  should  be  noted  that  an  error  is  introduced  in  this  approximation 
procedure  by  the  fact  that  the  differential  strips  of  cylinder  surface  overlap  in  the  vicinity  of  the 
poles  of  the  sphere  as  the  cylinder  is  rotated.  This  overlap  is  concentrated  in  the  “polar”  regions 
of  the  sphere,  however,  and  if  the  radius  of  the  sphere  is  large  the  NU  currents  will  decay  rapidly 
away  from  the  shadow  boundary  so  that  the  error  attributable  to  the  overlap  of  the  differential 
strips  can  be  expected  to  be  small. 

In  Figure  20  we  show  the  scattered  H-plane  far  field  for  a  sphere  of  size  ka  =  60  illuminated 
by  a  plane  wave  obtained  by  the  method  of  ILDC’s,  along  with  the  exact  field  obtained  from 
the  Mie  series.  The  back  and  forward  scatter  directions  correspond  to  0  =  0°  and  6  =  180°, 
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Figure  20:  Amplitude  of  scattered  H-plane  far  field  of  sphere,  ka  =  60,  illuminated  by  a  plane 
wave;  0°  <  0  <  180°. 

respectively.  The  detail  of  Figure  20  for  130°  <6<  180°  is  shown  in  Figure  21.  It  can  be  seen 
that  a  significant  improvement  upon  the  accuracy  of  the  PO  scattered  field  is  obtained  by  adding 
the  integral  of  the  circular  cylinder  NU  current  ILDC  to  the  PO  field  of  the  sphere.  A  similar 
result  was  obtained  for  the  E-plane  scattered  far  field  pattern.  Calculations  for  spheres  smaller 
than  ka  =  60  have  shown  that  the  ILDC’s  significantly  enhance  the  accuracy  of  the  PO  far  field 
for  spheres  as  small  as  ka  =  20. 


5  Conclusions 

The  nonuniform  currents  excited  in  the  vicinity  of  the  shadow  boundary  of  a  PEC  object  can 
contribute  significantly  to  the  total  scattered  far  field  of  the  object.  The  method  of  ILDC’s  is 
an  important  technique  for  calculating  this  NU  far  field  which,  added  to  the  PO  field,  gives  a 
considerably  more  accurate  approximation  to  the  total  scattered  far  field  than  does  the  PO  far 
field  alone.  The  approach  taken  in  this  report  to  obtain  ILDC’s  for  the  shadow  boundary  of 
a  PEC  2D  convex  cylinder  of  general  cross-section  is  to  integrate  the  product  of  the  free-space 
Green’s  function  and  the  NU  currents  over  a  strip  of  surface  of  the  2D  cylinder  of  differential 
width,  transverse  to  the  shadow  boundary.  Accurate  and  rapidly  calculable  Fock  approximations 
to  the  NU  currents  are  obtained  and  used.  For  circular  and  parabolic  cylinders  these  NU  current 
approximations  are  validated  numerically.  The  integration  required  to  obtain  the  ILDC’s  is  itself 
performed  rapidly  by  using  quadratic  polynomial  approximations  of  both  the  amplitude  and 
unwrapped  phase  of  the  integrand  of  the  radiation  integral,  thereby  enabling  the  integration  to  be 
performed  in  closed  form  involving  Fresnel  functions.  Examples  are  given  of  the  approximations 
of  both  the  NU  currents  and  the  integrand  of  the  radiation  integral.  As  an  application  of  this 
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Figure  21:  Amplitude  of  Scattered  H-plane  Far  Field  of  Sphere,  ka  —  60,  Illuminated  by  a  Plane 
Wave;  130°  <0<  180°. 

method  of  obtaining  ILDC’s  for  convex  PEC  cylinders,  the  ILDC’s  of  a  PEC  circular  cyhnder  are 
obtained  and  integrated  over  the  shadow  boundary  of  a  PEC  sphere  illuminated  by  a  plane  wave 
to  obtain  the  NU  current  far  field.  This  field,  added  to  the  PO  far  field,  is  shown  to  considerably 
ijjiprove  upon  the  accuracy  of  the  PO  approximation  of  the  exact  total  scattered  far  field  of  the 
sphere. 
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Appendix 

Closed-Form  Expressions  for  an  Integral 

In  this  appendix  we  present  closed-form  expressions  for  the  integral 

/(co,Ci,C2,C3,C4,C5,a,6)  =  /  (co5^ -I- Ci5 -f  (70) 

Ja 

occurring  in  the  expression  (61).  We  assume  that  cq,  Cj,  and  are  complex  and  that  C3,  C4,  C5, 
a,  and  b  are  real. 

Case  I:  C3  >  0 

Integration  by  parts  shows  that 
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where 


and 
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(73) 

(74) 


is  a  Fresnel  function.  When  |xo|  and  |xj,|  are  large  one  can  substitute  the  large- argument  ap¬ 
proximation  for  the  Fresnel  function  [20,  pp.  301-302] 
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into  (72)  to  get 


A  = 


i 

y/2'KC3 


— ) 

TT^xl) 


g— 4(036^  4-04  (>+cs) 


J_  _J _ l_\ 

Xa  TTXI  'k'^xI) 


g-z(c3a2+c4a+C5) 


(76) 


Numerical  calculations  show  that  it  is  advantageous  to  use  (76)  when  |x„l  >  10  and  jx^]  >  10. 
For  |xa|  <  10  or  jx^]  <  10,  the  routine  “frenel”  of  Numerical  Recipes  [21,  pp.  248-250]  can  be 
used. 
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Case  II:  C3  <  0 


The  integral  can  in  this  case  be  obtained  from 

/(Co,Ci,C2,C3,C4,C5,C,fe)  =  [/(Cq,  C*,  Cj, -C3, -C4, -Cg,  c,  6)]*  (77) 

where  =i=  indicates  complex  conjugation.  The  integral  on  the  right  side  of  (77)  can  be  computed 
by  the  formulas  given  in  Case  I  above. 

Case  III:  0  <  [csl  <C  1 

In  this  case  it  is  advantageous  to  obtain  I  from 

dl 

?  ^2 5  ^3?  ^4?  ^5^  5  ^2?  ^4?  ^55  ^5  “t”  (^0?  ?  ^25  ^ 7  ^4?  ^5 ? 

OCs 

where  ^ 

^(co,  Cl,  C2, 0,  C4,  C5,  a,  b)  -  -i  I  s^{coS^  +  ciS  +  C2)  (79) 

UC^  Ja 

The  integral  can  be  easily  evaluated  in  closed  form  by  integrating  by  parts  or,  more  simply,  by 
using  a  symbolic-integration  computer  program. 

Case  rV:  C3  —  0,  C4  0 

A  straightforward  calculation  shows  that 
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(80) 


Ccise  V:  C3  =  0,  C4  =  0 

A  straightforward  calculation  shows  that 

/(co,ci,C2,0,0,C5,a,6)  =  ^co(6^  -  a^)  +  ^Ci(6^  -  a^)  +  02(6  -  a)  e"*'".  (81) 

Equation  (66)  implies  that  the  quantities  and  e“*(=36^+c4i>+cs)^  occurring  in  many 

of  the  expressions  of  this  section,  simply  equal  and  respectively. 
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MISSION 

OF 

ROME  LABORATORY 


Mission.  The  mission  of  Rome  Laboratory  is  to  advance  the  science  and 
technologies  of  command,  control,  communications  and  intelligence  and  to 
transition  them  into  systems  to  meet  customer  needs.  To  achieve  this, 
Rome  Lab; 


a.  Conducts  vigorous  research,  development  and  test  programs  in  all 
applicable  technologies; 

b.  Transitions  technology  to  current  and  future  systems  to  improve 
operational  capability,  readiness,  and  supportability; 

c.  Provides  a  full  range  of  technical  support  to  Air  Force  Materiel 
Command  product  centers  and  other  Air  Force  organizations; 

d.  Promotes  transfer  of  technology  to  the  private  sector; 

e.  Maintains  leading  edge  technological  expertise  in  the  areas  of 
surveillance,  communications,  command  and  control,  intelligence,  reliability 
science,  electro-magnetic  technology,  photonics,  signal  processing,  and 
computational  science. 


The  thrust  areas  of  technical  competence  include:  Surveillance, 
Communications,  Command  and  Control,  Intelligence,  Signal  Processing, 
Computer  Science  and  Technology,  Electromagnetic  Technology, 
Photonics  and  Reliability  Sciences. 


